Computational spectroscopy of helium— solvated molecules: 
effective inertia, from small He clusters toward the nano-droplet regime 
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Accurate computer simulations of the rotational dynamics of linear molecules solvated in He 
clusters indicate that the large-size (nano-droplet) regime is attained quickly for light rotors (HCN, 
CO) and slowly for heavy ones (OCS, N2O, CO2), thus challenging previously reported results. 
Those results spurred the view that the different behavior of light rotors with respect to heavy 
ones — including a smaller reduction of inertia upon solvation of the former — would result from the 
lack of adiabatic following of the He density upon molecular rotation. We have performed computer 
experiments in which the rotational dynamics of OCS and HCN molecules was simulated using a 
fictitious inertia appropriate to the other molecule. These experiments indicate that the approach 
to the nano-droplet regime, as well as the reduction of the molecular inertia upon solvation, is 
determined by the anistropy of the potential, more than by the molecular weight. Our findings are 
in agreement with recent infrared and/or microwave experimental data which, however, are not yet 
totally conclusive by themselves. 

PACS numbers: 36.40.-c, 34.20.+h, 67.40.Yv, 36.40.Mr, 02.70.Ss 



INTRODUCTION 

The recent discovery that small molecules solvated in 
He clusters display a sharp rotational spectrum with well 
resolved lines is a spectacular manifestation of He su- 
perfluidity in such extreme confined conditions 0, 
Recent developments in quantum many-body simulation 
techniques are allowing for the determination of the first 
few excited states of interacting boson systems 0, 0, @ , 
as well as for an understanding of the relations existing 
between structure, dynamics, and superfluidity in these 
systems @0. In particular, it has now become possible 
to predict the dependence of the effective molecular rota- 
tional constant, B, on the cluster size, N, for clusters up 
to a few tens of He atoms || . 

Results based on the so-called projection operator 
imaginary-time spectral evolution (POITSE, an imple- 
mentation of diffusion quantum Monte Carlo aimed at 
estimating excitation energies) provide the following sce- 
nario for the dependence of B upon N 

0SGi3 : i) the 

inertia of heavy rotors is reduced by a factor 2-3 upon sol- 
vation, whereas the reduction is of the order of a few tens 
percent only in the case of light molecules; ii) the inertia 
of solvated heavy rotors (such as, e.g., OCS or SFg) would 
reach the large-size (nano-droplet) regime well before the 
first solvation shell is completed, whereas the convergence 
for light rotors would be much slower. This scenario was 
partially rationalized using the concept of adiabatic fol- 
lowing. Heavy molecules rotate slowly when excited and 
they usually have a stronger and more asymmetric in- 
teraction potential with the surrounding He atoms. As 
a result of the combination of these two effects, a frac- 
tion of the He solvent density would be dragged along 
by the rotation of the molecule, thus contributing to in- 



crease its inertia. This effect would be greatly reduced in 
light rotors both because they rotate faster, and because 
their interaction with He atoms is usually more spheri- 
cal. Both these facts hinder the dragging of He atoms, in 
qualitative agreement with experiments which positively 
support the first of the above two findings. The fact that 
only a fraction of the He density in the first solvation 
shell would adiabatically follow the molecular rotation 
would explain the alleged fast convergence of the molec- 
ular inertia to the nano-droplet regime for heavy rotors 
0, 0| . The alleged slow convergence for light rotors 
would instead be due to the interaction between molec- 
ular rotation and bulk-like collective excitations that is 
supposed to occur for large enough clusters 0, E3. |2(| . 
Experimental results often exist only for large clusters. 
Some results on small clusters begin to be available, but 
data in the intermediate regime are not yet sufficient to 
validate this picture. 

In this paper this problem is re-addressed using Rep- 
tation Quantum Monte Carlo (RQMC), a method that 
we believe to be the most appropriate existing to date 
to study the dynamical properties of interacting boson 
systems 0, 0- As paradigmatic cases of heavy and 
light solvated molecular rotors we choose OCS@Hcat and 
HCNQHeAr, respectively. In order to disentangle the 
role of the bare molecular inertia from that of the po- 
tential anisotropy, we have performed computer experi- 
ments in which the rotational dynamics of OCS and HCN 
molecules was simulated using a fictitious inertia appro- 
priate to the other molecule. These experiments indicate 
that the approach to the nano-droplet regime, as well as 
the reduction of the molecular inertia upon solvation, is 
determined by the anistropy of the potential, more than 
by the molecular weight. 
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Our paper is organized as follows. In Sec. II we briefly 
review our theoretical framework. Sec. III. A contains our 
results for OCS@Hejv which considerably extend the size 
range considered in a previous publication 0. We also 
address the problem of the accuracy of the molecule-He 
interaction potential \Ti\ and we show that avoiding the 
morphing procedure followed to brush up the potential, 
allows for a much better (actually, almost perfect) agree- 
ment between simulations and rotational spectroscopy 
experiments. In Sec. III.B we present results of simi- 
lar simulations made for HCNQHcat. In Sec. III.C we 
present computer experiments done on molecules with 
fudged masses (i.e. by using the inertia of one molecule 
and the interaction potential of the other), and we dis- 
cuss the implications of these results on the main issue 
raised in the present paper (i.e. the relative importance 
of the bare molecular inertia and of the anisotropy of the 
molecule-atom potential in determining the rotational 
properties of solvated molecules, as a function of the clus- 
ter size). Finally, Sec. IV contains our conclusions. 
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THEORY 

We consider a realistic Hamiltonian in which N He 
atoms are treated as point particles interacting via a pair 
potential, and the dopant molecule (either OCS or HCN) 
as a rigid linear rotor, with only translational and ro- 
tational degrees of freedom. Both the He-He and the 
He-molecule interactions are parametrized after accurate 
ab-initio quantum chemistry calculations [lllll2T , ll3j . For 
the He-OCS potential energy surface (P ES) , in particu- 
lar, we use the unmorphed version of Ref. I 111 In the next 
section we show that this version — which results from a 
direct fit to quantum-chemistry calculations — is consid- 
erably more accurate for a wide range of cluster sizes 
H than the morphed one which was refined so as to im- 
prove the predictions for the spectrum of the OCS@Hei 
complex . Contour plots of the He-HCN and the He- 
OCS PES arc displayed in Fig.Q showing a much greater 
strength and anisotropy of the latter. 

We simulate the system using the RQMC method 
0, Q- Details of the calculations are similar to those 
of previous studies for OCS@Hejv and CO@Hejv clus- 
ters 0, IhI ] . Briefly, the imaginary-time evolution oper- 
ator projects a state defined by a trial function $o onto 
| $(/?)) = exp(— /37J)|$o), which in turn approaches the 
exact ground state, Vf^, as (5 goes to infinity. Expectation 
values on |$(/?)), ($(/?)|0|$(/3)}/($(/3)|$(/3)), arc com- 
puted using a discretized path-integral representation of 
the imaginary-time evolution, which becomes exact in 
the limit of zero time step. The resulting paths are sam- 
pled with an efficient generalized Metropolis algorithm. 
The average (O) of the values taken by the operator O 
on the sampled paths gives an estimate of the desired 
expectation value within a known statistical error and 



FIG. 1: Potential energy surfaces of He-OCS (up) and He- 
HCN (down) used in this work. The He-OCS PES is the 
unmorphed version from Ref. ^] (see text). Contour levels 
start from V — to negative values, and they are spaced by 
5 °K. The He-OCS interaction potential has a deep equatorial 
doughnut-shaped well around the C atom, and two secondary 
polar minima in correspondence of the O and S atoms, the 
one at the S pole being broader and deeper. The He-HCN 
interaction potential, instead, has a shallow main minimum in 
correspondence of the hydrogen molecular pole and an even 
shallower secondary one slightly below the equator on the 
nitrogen side, the latter being barely visible in the figure. 



without any mixed-estimate nor population-control bi- 
ases. Here O can be either a static operator (such as the 
Hamiltonian itself or the atomic number density distribu- 
tion) or a time-dependent one (such as those needed to 
calculate time-correlation functions: Aexp(— tH)B. A 
single path thus entails the imaginary-time evolution of 
the system for a total time 2/3 + r, which in our sim- 
ulations is typically between 1 and 2 inverse K, with a 
discretization time step, e, of 10~ 3 K _1 . By explicitly 
checking the convergence of the results to the limits of 
large projection time and zero time step, we make sure 
that the computed quantities accurately represent their 
ground state values. 

The trial function is taken of the Jastrow form, 



$o = exp 



JY 



N 



i<j 



(1) 



where is the position of the i-th. atom with respect to 
the center of mass of the molecule, r, = |r,-|, 0, is the 
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angle between the molecular axis and r,, and is the 
distance between the i-th. and the j-th He atoms. On 
account of the anisotropy of the PES, the He- molecule 
pseudo-potential, U\, is expressed as a sum of products of 
Legendre polynomials times radial functions. Retaining 
five or six terms in the sum is enough to represent even 
the most anisotropic situation, namely a ring of up to 
5 He atoms tightly bound around the OCS molecular 
axis in correspondence of the principal minimum of the 
PES. The choice of the radial functions is not critical, 
as soon as their variational parameters ensure sufficient 
flexibility in the relevant range of Hc-molccule distances. 
The resulting variational parameters, whose number is 
in the range of several tens, are optimized by minimizing 
the variational energy with the variational Monte Carlo 
method. 

Information on the rotational excitations of the cluster 
can be obtained from the imaginary time autocorrelation 
functions 

c J (r) = (P J (n(r).n(0))) J (2) 

where n is the molecular orientation versor, and Pj is the 
Legendre polynomial of degree J. A generic imaginary- 
time correlation function can be represented as a linear 
combination of exponentials, whose decay constants are 
the excitation energies of the system, and whose coef- 
ficients are the corresponding oscillator strengths. The 
usefulness of this representation for obtaining dynamical 
properties from quantum simulations is in general rather 
limited. This is so because the calculation of a spectrum 
from imaginary-time correlations entails carrying out an 
inverse Lajjlace transform, a notoriously ill-conditioned 
problem In the present case, however, the situa- 

tion is not as bad because very few excited states con- 
tribute to cj{t). In fact, if the solvated molecule were 
isolated, only one rotational state would contribute to 
c,/(t) which would therefore have the form of a single 
exponential. Furthermore, the interaction between the 
solvent matrix and the soluted molecule is rather weak, 
so that this single-exponential picture is only slightly per- 
turbed. Last, and most important, the bosonic nature of 
the quantum solvent determines a low density of low- 
energy excitations. As a consequence of the scarcity of 
low-lying excitations available to couple with the molec- 
ular rotation, it is thus possible to reliably extract their 
energies by a multi-exponential fit. A more general analy- 
sis, based on the maximum entropy method, gives equiv- 
alent results 

RESULTS AND DISCUSSION 

Carbonyl sulfide 

Carbonyl sulfide (OCS) appears to have a special 
role in the spectroscopy of He-solvated molecules. OCS 



was used as a microscopic rotor to probe the superflu- 
idity of He nano-droplets in a celebrated microscopic 
Andronikashvili experiment 0. As a consequence of 
the paucity of low- lying excitations in the bosonic host, 
molecular rotational lines are exceedingly sharp and well 
resolved when the molecule is solvated in 4 Hc, while they 
are broad and unresolved in 3 He. In the former case the 
rotational spectrum is well represented by a quasi-rigid 
free-rotor model: 

Ej = BJ(J +1) - DJ 2 (J +1) 2 , (3) 

where B is the effective rotational constant of the 
molecule — renormalized by the interaction with the host 
to a value smaller but of the same order as in the gas 
phase — and D is the rotational distortion constant which 
is several orders of magnitudes larger than for the free 
molecule. 

A few years later, size-resolved rotational and vibra- 
tional spectra of small (N = 1 -r 8) OCS@HeAr clusters 
appeared 0] , indicating that fingerprints of superfluidity 
can be identified in clusters as small as TV = 8. These 
findings spurred a burst of theoretical research aimed at 
elucidating the relations existing among structure, dy- 
namics, and superfluidity in these small clusters jfj, l]]- 
A semi-quantitative agreement between POITSE results 
from diffusion quantum Monte Carlo simulations and 
experiments was reported in Ref. in the size range 
N = 1 -j- 8 (a minimum in the rotational constant was 
predicted for N = 6, while experiments show a decrease 
up to the largest resolved size, N = 8), while for N > 6 a 
fast convergence to the nano-droplet limit was observed 
(full convergence was claimed to be achieved for N = 20). 
A better agreement was found in Ref. (i. in that the 
minimum of the rotational constant (located at a cluster 
size N = 8) was found to be compatible with spectro- 
scopic data. Results were also presented for the rota- 
tional distortion constant, D, which was found to reach 
a minimum for N = 5, in agreement with experimen- 
tal findings |2j]. In addition, a detailed analysis — based 
on a comparison between the structure of the clusters 
and the He-OCS angular-current correlations — allowed 
to shed light on the microscopic mechanisms responsible 
for superfluidity. 

Appraising the quality of intermolecular potentials 

Much of the discrepancy between the results of quan- 
tum Monte Carlo simulations and experiments (as well 
as among different simulations, we should add) may be 
due to the quality of the He-molecule potential utilized 
for the simulations [|| . In Ref. ^] it was in fact shown 
that the difference between the predictions of Refs. and 
is indeed (almost) entirely due to the poorer quality of 
the potential utilized in Ref. [3. In Fig. H we com- 
pare the atomic density distributions of OCS@Hes and 
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FIG. 2: Upper panel: potential energy surface of He-OCS 
(we display here the morphed version from Ref. which was 
used in our previous study Q)- Lower panels: contour plots 
of the helium density profiles of OCSOHeg and OCS@Hei . 
For the density profiles, contour levels start from 0.001 with 
increments of 0.005, in units of A~ J . 

OCS@Hci with the OCS-Hc potential used in Ref. H 
According to the analysis of Ref. 0, the minimum in the 
effective inertia of the solvated OCS molecule occurs at 
the largest cluster size at which quantum tunneling be- 
tween the main, equatorial, and the two secondary, polar, 
potential wells is hindered by the energy barriers which 
separate them. It is clear that the larger the barrier, the 
smaller the tunneling, and the larger the corresponding 
effective inertia will be. In Ref. |6| it was indeed found 
that fudging the OCS-He PES — so as to enhance the 
potential energy barrier which separates the main well 
from the molecular poles — hardly affects the rotational 
constant for N < 5, while it increases the inertia for 
N = 6,7, 8, thus bringing the results of the simulations 
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FIG. 3: Effective rotational constant of OCS@Hejv, as a 
function of the cluster size, N. Full diamonds: results from 
RQMC simulations performed using the unmorphed potential 
of Ref. (see text). Open circles: experimental data from 
Ref. 0- The horizontal dashed line indicates the nano-droplet 
limit measured in Ref. 0- 

in much better agreement with experiments. It was later 
found that this too low a value for the relevant energy 
barriers is in fact an artifact of the morphing procedure 
adopted in Ref. 0, so as to bring an already very ac- 
curate OCS-He potential obtained from coupled-cluster 
quantum-chemical calculations into an even better agree- 
ment with the known spectra of the OCS@Hei com- 
plex. Such spectra depend on the details of the poten- 
tial around the minimum more than on the height of the 
barriers which separate incquivalcnt minima. The latter, 
instead, are expected to crucially affect the properties of 
the cluster when the main potentials well starts to be 
filled by He atoms (JV > 5, in the present case). We con- 
clude that the spectra of clusters with many He atoms — 
which depend on the PES far from the minimum — are a 
much more sensible benchmark of the quality of the PES 
than those of the helium-molecule dimer. A similar con- 
clusion was drawn in Ref. ^] where a different He-OCS 
PES was proposed, such that the rotational constants 
calculated from POITSE simulations agreed well with ex- 
perimental data available for N < 8. RQMC simulations 
performed with this PES in the size range 10 < N < 30 
show a behavior of the rotational constant very similar 
to that reported in Fig. |3 obtained from the unmorphed 
version of the PES of Ref. EH 

although the resulting val- 
ues are slightly smaller for the largest sizes (hence closer 
to the experimental nano-droplet regime: see below). 

Effective rotational constants 

In Fig. |31 we report the dependence of the effective ro- 
tational constant, B, on the cluster size, N, as calculated 
with the unmorphed He-OCS potential of Ref. 0. The 



5 



nano-droplct limit resulting from measurements on clus- 
ters of w 6000 atoms is indicated by a dashed line. 
Our results compare favorably with experimental data 
which are available up to N = 8 2J. In the small-size 
regime (N < 8) the rotational constant decreases with 
increasing cluster size because the He atoms trapped in 
the (equatorial) and in the first polar (near oxygen) well 
are dragged along by the molecular rotation, thus in- 
creasing the effective inertia. For N > 8 — as the second 
polar (near sulphur) well starts to be filled — an increasing 
fraction of the He atoms can freely tunnel among differ- 
ent wells, thus not contributing to the molecular inertia. 
The quantum fluid nature of the He solvent is such that 
tunneling is a collective process in which not only the 
excess atoms not fitting in the main potential well take 
part, but even those which are tightly bound to it. For 
this reason, once tunneling among different wells is made 
possible by the spill-out of excess atoms, this process de- 
termines a decrease of the molecular inertia, i.e. an in- 
crease of the rotational constant for an increasing cluster 
size. As a matter of fact, it was shown in Ref. that 
the value of the correlation between the molecular angu- 
lar momentum and the atomic angular current — which is 
maximum in the main equatorial potential well — starts 
decreasing as the rotational constant increases past the 
minimum (N > 8). This decrease continues until the 
first solvation shell is completed at a cluster size N pa 20, 
around which the rotational constant seems to stabilize. 
As the second solvation shell starts to build, however, 
quantum exchange cycles involving atoms from this shell 
would contribute to further, although weakly, decrease 
the molecular inertia. A similar behavior of the evolution 
of the rotational constant for sizes shortly beyond com- 
pletion of the first shell is also observed in simulations of 
clusters doped with N2O and CO2 (lflllflj. which are all 
molecules having a qualitatively similar interaction with 
He atoms. Our findings demonstrate that — contrary to 
a commonly accepted assumption — He atoms from outer 
(larger than the first) solvation shells do affect the molec- 
ular inertia. For N2O and CO? this is also supported by 
experimental evidence |2fl l2l| : the measured value of 
B for the largest cluster with secure assignment of spec- 
tral lines (N = 12 for N 2 and N = 17 for C0 2 ) is 
significantly higher than the nano-droplet limit, with no 
plausible signs of convergence within the first shell. A 
role of outer shells in bringing down the effective rota- 
tional constant is thus to be expected. How far from 
the molecule does this effect extend, our simulations — 
which arc limited at present to a few tens atoms — cannot 
say yet. This finding is at least compatible with current 
phenomenological models of the inertia of He-solvated 
molecules, which predict a lower contribution of the sol- 
vent to the molecular inertia, with decreasing atomic den- 
sity. In hydro-dynamical models a lower atomic density 
would determine a reduced kinetic energy of the irrota- 
tional flow of the solvent; in a two-fluid model, instead, a 
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FIG. 4: Radial density profile p(r) of He in OCSOHejv for 
TV = 20 (dashed line) and N = 50 (solid line). 

lower inertia would simply follow from a reduction of the 
non-superfluid component of the atomic density. Com- 
ing down from the nano-droplct regime, a decrease of the 
atomic density in the first solvation shell(s) is indeed ex- 
pected, as a consequence of the reduced pressure exerted 
by the outer shells. This is clearly demonstrated in fig. 01 
in which the He radial density profiles around OCS are 
compared for N = 20 and N = 50. The radial den- 
sity in the first shell is significantly lower for the smaller 
cluster. On purely classical grounds, a competing effect 
could arise if the spatial extension of the first shell was 
sensitive to the pressure release, thus affecting the sec- 
ond moments of the atomic density. The density profiles 
shown in Fig. 0] however suggest that this effect is small, 
since the positions of the peaks in the first shell hardly 
change between clusters of 20 and 50 He atoms. 

Our results demonstrate that in OCS@Hejv clusters 
the effective rotational constant does not attain its 
asymptotic limit upon completion of the first solvation 
shell, being in fact higher at this size. When the second 
shell starts to build up, the value of B further increases 
(arguably, via quantum exchange cycles involving parti- 
cles of both shells 0,H2). Whether, in larger clusters, 
the missing inertia will be recovered by a change in the 
density of the first shell or by a direct contribution from 
the outer shells, or both, remains to be investigated. 



Hydrogen cyanide 

Hydrogen cyanide (HCN) is considered as a prototype 
of light helium-solvated rotors 0, 0| , for which the as- 
sumption of adiabatic following breaks down. This has 
been demonstrated both experimentally, by comparing 
the rotational constants of HCN and DCN in the nano- 
droplet regime |23l |. and theoretically, by comparing the 



FIG. 5: Contour plot of the incremental density, Apjv(r) = 
pjv(r) - p N -i(r), for selected HCN@Hejv clusters. The HCN 
molecule has its center of mass at the origin, with the hydro- 
gen atom on the positive x direction. Distances are in A, and 
contour levels start from 0.0001 with increments of 0.001, in 
units of inverse cubic A. 
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FIG. 6: Chemical potential, fi(N) = E(N) - E(N - 1), for 
HCNQHejv clusters. The dashed line is the bulk limit. 



He density profiles obtained from simulations performed 
including and, in turn, neglecting the molecular degrees 
of freedom |24j. The calculation of the molecular rota- 
tional constant, B, in large droplets thus defies the appli- 
cation of either the two-fluid |9[ and hydrodynamic mod- 
els |25j . which both rely on adiabatic following, albeit 
in a different manner. On the other hand, direct calcu- 
lation of rotational excitations by the POITSE method 
jiJl indicated that the nano-droplet value would not be 
reached even for N = 25, well beyond completion of the 
first solvation shell. This slow convergence of the ro- 
tational constant as a function of the cluster size was 
later attributed to the coupling of molecular rotation to 
phonon-like excitations of the solvent which would de- 
velop only in the nano-droplet regime |26l | . Unfortunately 
experimental data for the effective rotational constant of 
HCN are only available, so far, for very large droplets 
|2fl| . We will see however how a comparison of our the- 
oretical results for HCHQHcjv with both theoretical and 
experimental results which are available for the closely 
related CO@Hejy system will allow to draw a number of 
important and non trivial conclusions. 



Structural properties 

In Fig. [J]we report the incremental atomic density pro- 
file of HCN, Apjv(r) = pn{^) — /?jv_i(r), as calculated 
for a few selected cluster sizes. For the binary complex, 
HCN@Hei , the He density forms a broad cap around the 
H end of the HCN molecule, centered on the main well 
of the PES. As N increases, the He density smoothly 
piles up with no visible signature of the secondary, sub- 
equatorial, minimum. The incremental density is mainly 
localized on the H side of the HCN molecule for N up to 6, 
mainly on the N side for N between 8 and 14, and nearly 



isotropically thereafter. Starting from N around 16, the 
incremental density shifts toward larger distances from 
the HCN center of mass (see the inner and the outer con- 
tour levels in Fig.[SJ). We consider this behavior as an in- 
dication that the first solvation shell is completed around 
N w 15, although shell effects are by no means sharp in 
this system (for instance, much more pronounced effects 
are seen at N = 12 in clusters of para- hydrogen seeded 
with HCN Hill). 

The dependence of the cluster ground-state energy, 
E(N), on the size of the system reflects the extremely 
smooth evolution of the He density profile, as shown 
in Fig. [fj] which reports the chemical potential, /i( N) = 
E(N) - E(N - 1). For TV up to a dozen the effects of the 
He-HCN interaction dominates: the chemical potential 
stays almost constant, decreasing only very weakly with 
the cluster size, as a consequence of the He-He interac- 
tion. For 12 < TV < 20 the first solvation shell is filled 
and the kinetic energy of He atoms increases due to their 
closer packing, thus determining a rise of the chemical 
potential. For TV > 20 fi(N) is stabilized again at value 
which is higher than the bulk limit, due to the smaller 
effects of the He-He interaction in this size range. 

The fact that the first solvation shell is completed in 
the size range 12 < TV < 20 is confirmed by an analysis of 
the angular correlations of He atom pairs. Let us define 
the He-He pair distribution of the dihedral angle with 
respect to the molecular axis, <j>, as: 

^ = Nl^(E S ^-^ + ^)- ( 4 ) 

C(4>) gives a measure of the (weak) tendency of the He 
atoms to cluster together on the same side of the molecu- 
lar axis, due to the He-He attraction 0|. The integrated 
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FIG. 7: Integrated pair distribution of the dihedral angle, de- 
fined in Eq.|K| as a function of the cluster size, for HCNQHejv. 



quantity, 



r/2 

M = / C(0)d0 - 1/2, 
/o 



(5) 



displayed in Fig. [7| as a function of N, reveals that this 
tendency is smallest around N = 15, corresponding to 
completion of the first solvation shell where, for steric 
reasons, the angular distribution of He atoms around the 
molecule is most uniform. 



Rotational excitations 

All the above features of He-solvated HCN molecules 
are similar to the CO case, already studied in a previous 
work 0| . In Fig. [5]we report the energies of the two low- 
est rotational excitations with angular momentum J = 1 , 
along with the corresponding spectral weights. We note 
that our predictions for N > 10 considerably differ from 
those of Ref. which were obtained with the POITSE 
method. We cannot offer any explanation for this dis- 
crepancy. We can only observe that the value of B cal- 
culated in Ref. |^| for TV = 25 is even larger than in 
the gas phase, a fact that can hardly be explained on 
physical grounds. In fact, while the reduction of one of 
the moments of inertia below its value in gas phase — 
experimentally observed in the binary complex — can be 
explained by the considerable relative mass redistribu- 
tion occurring upon formation of the binary complex, 
similar effects are hard to justify at cluster sizes where 
the density profile evolves with N in a smooth and al- 
most isotropical way. This seems to suggests that the re- 
sults of Ref. may be affected by some inaccuracies for 
the largest sizes considered in that work. On the other 
hand, while it is difficult, for large cluster sizes, to en- 
sure full convergence with respect to projection time, as 
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FIG. 8: Upper panel: Rotational energies of HCNQHejv as 
functions of the cluster size, computed by RQMC (a-type, 
diamonds and b-type, filled circles); results from the POITSE 
calculations of Ref. |24| are indicated by squares. The dashed 
line is the effective rotational constant of HCN in the nano- 
droplet limit |23l|. while the empty circle at TV = is the gas- 
phase value. Energy units are cm -1 . Lower panel: Spectral 
weight of the a-type line (diamonds) and the fc-type line (filled 
circles) . 



well as full ergodicity in path sampling, we note that the 
present calculations for HCN arc indirectly supported by 
the agreement between the experiment |27j and similar 
calculations for the closely related CO-He clusters [Tij . 

Two series of excitations, called a-type and &-type 
lines, evolve smoothly from the known end- over- end and 
the free molecule rotational modes of the binary complex, 
respectively. The &-type line starts off with a stronger 
spectral weight, then it quickly weakens and eventually 
disappears. Note how the intensity of the &-type line fol- 
lows the decline of the He-density anisotropy displayed 
in Fig. The existence of two relevant spectral lines 
(instead of a single line appropriate to the spectrum of 
a linear rotor) is in fact a manifestation of the dynami- 
cal anisotropy of the atomic-density distribution around 
the molecular axis, which lowers the cylindrical symme- 
try of the molecular rotor. For N > 15 the rise of the 
He-density anisotropy does not give rise to any significant 
line splitting because it is due to atoms in the second sol- 
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FIG. 9: Comparison between the positions of the spectral 
lines reported in Fig. [SJwith the lowest mode in the spectral 
resolution of the correlation j(t). 



vation shell, which are only very weakly coupled to the 
molecular rotation. 

The coupling between the molecular rotation and He- 
density fluctuations is best understood using the (imagi- 
nary) time correlations of the versor, u, pointing from the 
molecular center of mass towards the center of mass of the 
complex of He atoms. Such a time correlation function, 
j(t) = (u(0) • u(t)), contains information on the ener- 
gies and spectral weights of the J = 1 cluster excitations 
whose character is predominantly that of a He-density 
fluctuation. In Fig. [§]we display the lowest-lying excita- 
tion energy extracted from the spectral resolution of "/(t) 
as a function of the cluster size, and compare it with the 
positions of the a-type and b-type lines already reported in 
Fig. |SJ We see that this density-fluctuation excitation — 
which has a relatively strong spectral weight — is degen- 
erate with the a-type line for N < 10. This mode has 
been interpreted as a cluster excitation in which part of 
the He density is dragged along by the molecular rotation 
|l6j . This interpretation is trivial for the binary complex 
(it corresponds to the end-over-end rotation), but it has 
no obvious visualization for larger clusters. For N > 10, 
the energy of this mode departs from the a-type line, 
and the b- type line correspondingly disappears. This re- 
lates to the onset of a situation where the He density is 
decoupled from the molecular rotation. 

We come now to the main concern of this paper, i.e. 
the convergence of the effective rotational constant to its 
asymptotic, large-size, value. Our results indicate that 
the residual renormalization of the effective inertia does 
not change significantly upon further growth of the clus- 
ter beyond, say, N = 15. This fact is in complete analogy 
with the findings reported in Rcf. for clusters doped 
with CO (and in agreement high resolution IR spectra 
recently obtained for He clusters seeded with CO up to 
N = 20 [13), but at variance with the results reported 



for HCN in Ref. |2J and with a commonly accepted 
view. With a similar rotational constant and a simi- 
lar interaction with He, CO and HCN are not expected 
to behave very differently upon solvation in He clusters. 
The experimental results for CO do not seem compat- 
ible with a large variation of B between, say, N = 15 
and 25, thus challenging the idea of slow convergence to 
the nano-droplet limit as a general feature of the effec- 
tive inertia of quantum solvated light rotors. Although 
a conclusive answer will require the measurement and 
assignment of spectral lines for even larger clusters, we 
believe that the agreement between our previous calcula- 
tions and high-resolution IR measurements for the closely 
related CO@Hejv system warrants a considerable trust in 
the present results and in the conclusions on the approach 
to the nano-droplet regime based on them. 

In conclusion, the results reported here for HCN and in 
Refs.Hilfor CO (the latter being supported by the experi- 
mental study of Ref. suggest that for these light rotors 
the asymptotic value of the effective rotational constant 
is reached well before completion of the first solvation 
shell. We will see in the next Section that this behavior 
is better attributed to the weak anisotropy of the poten- 
tial, rather than to the small value of B in gas phase. 



Fudged molecules 

In this section we consider the fractional reduction of 
the gas-phase rotational constant Bo upon solvation in 
He nano-droplets, A = B/Bq. The observed general 
trend (see e.g. Fig. 13 in Ref. Hoi) is that lighter ro- 
tors tend to have larger values of A. A suitably defined 
amount of adiabatic following \?>l\ has been proposed as 
the key physical property responsible for the value of A. 
Qualitatively, the analysis of Ref. supports the simple 
picture that both a small molecular inertia and a weakly 
anisotropic interaction lead to a large value for A. How- 
ever, there are several exceptions: for instance, N2O has 
Bo = 0.30 K and A = 0.17, whereas OCS (a typical en- 
try in the mentioned general trend) has B$ — 0.15 K and 
A = 0.36. Given that for OCS there is saturation to al- 
most full adiabatic following, |3l| the significantly smaller 
A value for N2O cannot be explained by resorting to the 
concept of adiabatic following alone. Therefore, it seems 
useful to gain further insight by disentangling the role 
of the PES anisotropy from that of the gas-phase in- 
ertia. To this purpose, we have performed simulations 
with two fictitious molecules: f-OCS (fudged OCS), with 
the PES as OCS and the same Bq value as HCN, and 
f-HCN (fudged HCN), featuring the HCN-He interaction 
and the Bq value of OCS. Note that the gas-phase rota- 
tional constants of OCS and HCN are in a ratio of about 
1:7. 

The rotational energies of f-OCS<Q?He7v are shown in 
Fig- El an d the corresponding spectral weights in Fig. 1111 
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FIG. 10: Rotational energies of f-OCS@Hejv as a function of 
the cluster size (a-type, filled circles, and fc-type, diamonds). 
The cross at TV = shows the fictitious value of 2Bo- Inset: 
the detail of the a-type line. 
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FIG. 11: Spectral weight of the a-type line (filled circles) and 
the fe-type line (diamonds) for f-OCS@Hejv. 



The general appcarcncc of both quantities for f-OCS is 
closer to HCN than OCS, due to the presence of an a- 
type and a &-type line, with the spectral weight of the 
latter decaying with increasing TV. However, at variance 
with HCN, the a-type line for large TV approaches an en- 
ergy significantly smaller than the gas-phase value, 2Bq. 
Taking the value at TV = 30 as an estimate (very likely an 
overestimate, see Section III. A) of the asymptotic value, 
we obtain A ~ 0.33, which is close to — and somewhat 
smaller than — the value of OCS. 

The a-type line of f-OCS shows a minimum at TV = 6 
or 7, i.e. at a smaller size than OCS. As shown in Fig. 1121 
already for TV = 8 there is a significant He density all 
around the molecule, which is presumably responsible 
for the turnaround of the B value (although the relation 
between density profiles and turnaround of the B value 
could be not so straightforward for f-OCS, due to the 
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Contour plot of the helium density profile of f- 



rcsidual spectral weight in the &-type line). Regardless 
of the implications for the turnaround, we stress that 
the significant difference between the density profiles be- 
tween f-OCS and OCS (see Fig.dand Fig.Hfor TV = 8) 
implies a significantly different amount of adiabatic fol- 
lowing. The fact that the A values of f-OCS and OCS are 
nevertheless very similar, indicates that the anisotropy of 
the potential, rather than the dynamical regime implied 
by the gas-phase inertia, is mainly responsible for the 
renormalization of the rotational constant upon solva- 
tion. 

Similar conclusions hold for f-HCN as well. The energy 
of the lowest rotational excitation with J = 1 is shown 
in Fig. EH as a function of TV. Already for TV = 3 the 
spectral weight of this excitation exceeds 90 percent, and 
in this respect f-HCN is closer to OCS than to HCN (i.e. 
the spectrum looks that of a linear rotor for TV > 3). 
However as far as the value of A is concerned, the effect 
of fudging the gas-phase inertia is very small. 

In order to estimate the asymptotic limit of the ro- 
tational constant, wc assume that it is given by the 
value at the largest-size cluster simulated (TV = 25), 
noting that the evolution of B(N) in Fig. ^5] is nearly 
flat for TV > 12 (for HCN, this assumption would give a 
very good agreement with the experimental nano-droplct 
value, see Fig. El . This gives A = 0.90 for f-HCN, close 
to — and somewhat higher than — the value 0.81 measured 
in HCN 0. 

Our findings indicate that the potential alone has a 
dominant role in establishing the value of A, at least 
in the range of physical parameters appropriate to the 
linear molecules studied here. Furthermore, the effect 
of reducing the molecular inertia while keeping the PES 
fixed is to slightly decrease the A value. This result, 
in agreement with an experimental study of HCN and 
DCN [23| . contradicts the empirically established corre- 
spondence between light rotational inertia and large A 
values. We conclude that such a correspondence is due, 
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FIG. 13: Rotational energies of f-HCN@Hejv as a function 
of the cluster size. The cross at N = shows the fictitious 
value of 2 Bo. 



to a large extent, to the generally small anisotropy of 
the interaction between light rotors and Helium an ac- 
cidental effect, as far as the rotational dynamics of the 
solvated molecule is concerned. Within this picture, the 
behavior of N 2 is not an anomalous case, but merely 
a consequence of the stronger stiffness and anisotropy of 
the N 2 0-He potential with respect to, say, OCS or C0 2 . 

We finally note that a previous calculation for SFg with 
a fictitiously small gas-phase inertia would support the 
opposite conclusion that, for given interaction with the 
solvent, a lighter molecule would have a larger value of 
A (3^ ■ While this result could shed some doubts on the 
generality of the conclusions drawn from the analysis of 
a few linear molecules, we believe that the calculations of 
Ref. probably deserve further analysis because they 
were obtained using a fixed-node approximation, whose 
accuracy is not warranted especially for light rotors [24| • 



CONCLUSIONS 

Quantum Monte Carlo simulations, in conjunction 
with high-quality inter-particle potentials, have reached 
a remarkable degree of accuracy for the calculation of the 
rotational dynamics of molecules solvated in He. Using 
the RQMC method, we have studied the evolution of the 
rotational excitations with the number N of solvent He 
atoms for a prototype heavy rotor, OCS, and a prototype 
light rotor, HCN (for OCS, preliminary results have al- 
ready been presented in Ref.Q). The size range explored, 
larger than presently attained with number-selective IR 
and/or MW spectroscopy, allows us to draw a series of 
conclusions on the approach of the rotational constant 
to its asymptotic value in the nano-droplet limit. Our 
results entail a substantial revision, both quantitative 
and qualitative, of the common view that the asymptotic 



limit would be essentially determined by the amount of 
adiabatic following and that — at least for heavy rotors — 
it would be reached well before completion of the first 
solvation shell. 

The rotational constant of OCS, after the undershoot 
of the nano-droplet value and the turnaround which sig- 
nals the onset of superfluidity, crosses again its asymp- 
totic limit at N = 12; moreover, starting with the be- 
ginning of the second solvation shell around N ps 20, 
it develops a further structure with an extremely broad 
maximum, and a possible hint of a (final?) decrease only 
seen at the largest size we studied, N = 50. This fea- 
ture, which parallels similar findings for C0 2 and N 2 0, 
definitely supports slow convergence to the nano-droplet 
value. A strikingly different behavior is found for HCN. 
In this case, a linear-rotor-likc spectrum is found for N 
larger than 10, and the resulting rotational constant stays 
constant in a wide range (say 15 to 50), with a value close 
to the measured value in the nano-droplet limit. This re- 
sult contradicts the expectation of a slow convergence to 
the asymptotic limit (determined by coupling of molec- 
ular rotation with well-developed bulk-like excitations of 
the solvent) as a general property of light rotors. In- 
deed, a behavior very similar to that illustrated here for 
HCN was found for another light rotor, namely CO. It 
would be tempting to propose that fast convergence to 
the nano-droplet limit is a general rule for light rotors. 
However, preliminary experimental results [33j seem to 
indicate that, for CO, the asymptotic limit is significantly 
lower than inferred from the nearly constant value of the 
rotational constant in the size range from a dozen to a 
few tens He atoms. The possibility of defining a general 
trend for light rotors thus deserves further investigation. 

In order to establish the relative importance of the bare 
molecular inertia and of the strength and anisotropy of 
the He-molecule interaction in determining the approach 
of the rotational dynamics to the nano-droplet regime, we 
have also performed computer experiments in which the 
molecular inertia was intentionally modified. To this end, 
we have considered two fudged molecular species, f-OCS 
and f-HCN, i.e. OCS and HCN with fictitiously small 
and fictitiously large values of the gas-phase rotational 
constant, respectively (appropriate in fact to the other 
molecule). Perhaps the most important feature which 
was attributed to the predominant role of the bare molec- 
ular inertia is the amount of renormalization of the gas- 
phase rotational constant, B, upon solvation: the strong 
(weak) reduction of B observed for heavy (light) rotors 
was attributed to the large (small) amount of adiabatic 
following. Our results indicate that the fractional reduc- 
tion of the gas-phase rotational constant upon solvation is 
slightly stronger for f-OCS than for true OCS, despite the 
obvious fact that adiabatic following is much larger for 
the latter. Likewise, the reduction calculated for f-HCN 
is somewhat weaker for f-HCN than for true HCN. The 
same trend was experimentally observed, with a smaller 
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variation of Bq, in a comparative study of HCN and 
DCN [2^]. This clearly shows that it is the strength and 
anisotropy of the He-molecule interaction, rather than 
the bare molecular inertia, which is mainly responsible 
for the renormalization of the rotational constant in the 
nano-droplet regime. In this perspective, the classifica- 
tion into heavy and light rotors thus retains its validity 
only to the extent that heavier molecules tend to have 
stronger and more anisotropic interactions with He. 

We thank Giacinto Scoles for many fruitful discussions. 
Part of the calculations presented in this work were done 
thanks to the Iniziativa Calcolo Parallelo of the Italian 
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